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ABSTRACT 

We revisit, via a very simplified set of equations, a linear streaming instability (tech- 
nically an overstabili ty), which is present in, and potentially important for, dusty 
protoplanetary disks ( Youdin and Goodman|[2005 1 . The goal is a better understand- 
ing of the physical origin of such instabilities, whicn are notoriously subtle. Rotational 
dynamics seem to be essential to this type of instability, which cannot be captured 
by one-dimensional Cartesian models. Dust 'pileups' in moving pressure maxima are 
an important triggering mechanism of the instability, and drag feedback of dust upon 
the gas allows these maxima to be strengthened. Coriolis forces and the background 
drift counteract the effects of the pressure force. 
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The formation of planets via dust accretion starting 
with micron-sized grains and ending with 10* — 10^-km-sized 
bodies in a protoplanetary disc is a multiphase process (see 
e.g. the review by Chiang and Youdin (20101, henceforth 



the Roche density. However, these mechanisms assume a low 
level of global turbulence. While the magnetorotational in- 



CYIO). Grains smaller than roughly a centimetre coagulate 
via sticking (e.g. Dominik et aZ.|2007 l, while for the largest 
sizes, gravity of the planetary bodies plays a dominant role 
(e.g. pa|[2010l ). 

Between these two limits, there is a gap that neither 
mechanism seems able to leapfrog. Beyond a size of 1 cm - 1 
m, collisions occur at velocities of order 1 — 10ms~^. These 



stability ( Balbus and Hawley 1998 1 may be suppressed in 
the so-called dead zone ( jGammiellO gB*), fluctuations arising 
from forced density waves or other instabilities could still 
be present ( Fleming and Stone|2003 Lesur and Papaloizou 



[2010 ) 



K** lead either to fragmentation or to bouncing (Zsom et al 
' i-H 

X 



rather than to systematic aggregation. Worse yet, 
particles of this size will drift sunward at the dramatic rate 
of 1 AU per century ( Weidenschilling|1977 1 . This is because 
solids tend to orbit the Sun at the full Keplerian speed, while 
the gas moves slightly slower because of the support from an 
outward directed pressure force. Thus, the solids experience 
a headwind that removes angular momentum via gas drag. 
This effect is the cause of the "metre-size barrier" . 

Self-gravity of the dust could in principle promote ac- 
cretion past this barrier provided some process enhances 
the dust density. A classic scenario for this is one in which 
the vertical component of the Sun's gravity causes sedi- 
mentation to the disc midplane and increases the dust den- 
sity ([Goldreich and Ward||1973| . Although [Weidenschilling| 
( 1980 1 had suggested that Kelvin-Helmholtz instabilities 



|2010[ [Latter et aZ.||2010| ). Concentrating the dust is prob- 
lematic. 

However, interesting effects already occur as soon as 
the dust-to-gas ratio is of order unity, say around the mid- 
plane, or as a result of turbulent concentration (e.g. [Cuzzi[ 
et aL|2003 l. This is because the dynamical back-reaction of 
the dust on the gas is then important, while self-gravity it- 
self may still be negligible. In fact, |Youdin and Goodman| 
( 2005 1 (hereafter YG05) discovered that the interpenetra- 



tion of dust and gas in a Keplerian disc was linearly over- 
stable: a process that they named the streaming instability, a 
convention we shall henceforth adopt. The dust particles can 
behave collectively via their interaction with the gas, lead- 
ing to disordered flows that facilitate grain clumping — as 



witnessed by the numerical simulations of [Youdin and Jo- 



hansen (20071 (hereafter YJ07) and Johansen and Youdin 
( [2007[ ). [Bai and Stontr ( 2010[ ) confirmed that in the nonlin- 
ear regime, dust densities up to three orders of magnitude 
above the background could be attained, which [Johansen] 
\et al.\ ( |2007| ) argued would suffice to instigate gravitational 
clustering. As relatively large (10 cm - 1 m) particles are 



could prevent the Roche density from being reached, re- 
cent studies indicate that this is not an issue for discs a 
few times more massive than the minimum mass solar neb- 
ula, or with supersolar metallicities (e.g. Chiang[[2008 Lee 
\et aZ.|2010'a|b| . Moreover, [Youdin] ( [20lT| showed that inclu- 
sion of gas drag could alleviate the need to actually reach 



needed for significant clumping ( Bai and Stone|2010 l, it re- 
mains to be demonstrated whether solids can have grown 
this large beforehand ( Lee et aZ.|2010b" l. 

Although much attention has been devoted to the 
streaming instability in the past years, its physical inter- 
pretation in the linear phase remains somewhat obscure. 
For example, as YJ07 point out, the idea that radial drift 
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slows down in overdense regions and leads to a sort of 'traffic 
jam', does not explain the onset of growth in the linear phase 
— such reasoning leads, in fact, to stable wave propagation. 
Such traffic jams, however, appear to be relevant for the non- 
linear phase (JY07). JY07 showed that drafting, analogous 
to pelotons in bicycle races, also was not necessary for the 
initial generation of particle overdensities, although it could 
enhance growth if present. CYIO discuss a one-dimensional 
vertically integrated toy model drawn from [Goodman and| 
Pindor (20001, in which a (single) fluid is subject to a "col- 
lective drag" acceleration in the sense that a drag coefficient 
depending upon the fluid density emerges from collective be- 
havior. This system is found to be overstable. It is unclear, 
however, in what sense the streaming instability conforms 
to this notion of collective drag since it is local in nature: 
the equations are not vertically integrated and the drag ac- 
celeration does not depend ab imtio on the density of the 
fluid acted upon. 

We are therefore motivated to reconsider the origins of 
dust clumping. In this paper, we revisit the linear stream- 
ing instability. In §2, we review the basic equations leading 
to the YG05 results and propose, in an appropriate limit, 
a markedly reduced system in which the instability arises. 
We make the point in §3 that rotational dynamics are es- 
sential to understand the streaming instability. In §4, we 
propose an interpretation of the streaming instability. In §5, 
we summarise our conclusions. 



2 STREAMING INSTABILITY AND 
REDUCED SYSTEMS 

In this section, by way of establishing notation (largely 
that of YG05), we briefly review the fundamental equations 
along with some of their more important results. We then 
propose a significant reduction of the system of equations 
that involves two successive levels of approximation, yet still 
produces the streaming instability. 



2.1 Fundamental two-fluid equations and local 
stability analysis 

Let us model the dust as a collection of identical, inde- 
structible spheres, of mass density pp and velocity Vp. The 
gas density and velocity, on the other hand, we denote by 
Pg and Vg, respectively. Their evolution equations are (see 
Appendix A for justification of the dust equations): 



dpp 
dt 

dpg 
dt 



+ V ■ (ppVp) = 



^stop 



V„Y„ = -!^^R + 



Pp Vp 



Ps 



VP 

Pa ' 



where SI is the Keplerian angular velocity. Dp 



(1) 
(2) 
(3) 

(4) 

d/dt + 



Vp,g ■ V the particle/gas Lagrangian derivatives, and R the 



cylindrical vector radius. We ignore the vertical component 
of the solar gravity, as well as self-gravity. 

The particle stopping time tstop depends on the size and 
velocity regime (e.g. ^Weidenschilling|1977 1. For example, for 
particles that are both small compared to the gas mean- 
free-path and drifting subsonically relative to the gas (as 
appropriate, e.g., for chondrule-sized bodies at a few AUs 



in a Minimum Mass Solar Nebula), Epstein's law (Epstein 



19241 applies: 



^stop 



PsO, 
PgVr' 



(5) 



with vt = \/9>kBT /-nm (with m the molecular mass and T 
the temperature), roughly the sound speed, ps the internal 
density of the grains (not to be confused with pp) and a the 
grain radius. Regardless of the relevant drag law, one de- 
fines a dimensionless stopping time measuring the coupling 
of dust to dynamical disturbances: 



(6) 



By virtue of Newton's third law, an all-important feedback 
term of the dust on the gas appears in the Euler equation for 
the gas (|4|. The system of equations is closed by assuming 
gas incompressibility (the Boussinesq approximation), as in 
YG05. 

It is convenient to express equations ([l])-(|4| in terms 
of centre-of-mass velocity V = (ppVp -I- Pg'Vg)/p (with p = 
pp + pg the total density) and the relative dust-to-gas drift 



AV; 



rg (YG05) 

dp 
dt 



+ V • (pV) = 



V • (V - ^AV) = 



dV 



+ V- VV + F(AV^ 



-n^R 



VP 

p 



(7) 



(8) 



(9) 



' V.V(AV)+(AV.V)V+G(AV^) = -J^^ + ^, 

Pg istop Pg 

(10) 



dt 

with: 



F(AV^ 



iv 



(11) 



G(AV^ 



Ps 



AV- V 



Ps 



AV 



Pp 



AV - V 



Pp 



AV 



(12) 

We shall now work in the so-called shearing-sheet approxi- 
mation: We neglect all curvature terms and disc-scale gradi- 
ents since the lengthscales of interest are much smaller than 
the heliocentric distance, or even the pressure scale height. 
We also adopt Cartesian coordinates in a frame corotating 
with the Keplerian fiow at a fixed heliocentric distance Po, 
at angular velocity Q. We denote hy x, y and z coordinates 
in the radial, azimuthal, and vertical (that is, perpendicular 
to the midplane) directions, respectively, and e^;. By and e^ 
the corresponding unit vectors. 



Using this model, [Nakagawa et al.\\ 1986 1 computed the 



equilibrium solution to equations 



■^(|1986 
7|-((l0f, 



v = (- 



20, 



(13) 
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AV = 



Qeistop 



+ 



fgSe^istop 

2(1 + 



(14) 



where = Pg.p/p the mass haction of gas and particles 
(of uniform densities in this equilibrium flow), and 

1 dP 
pdR 



9e 



(15) 



is the pressure-induced acceleration on the gas+dust fluid. It 
is the free energy stored in this gradient that the streaming 
instabihty can access (YG05). This gradient, consequently, 
introduces a lengthscale which helps determine the charac- 
teristic scale of the instability. We denote this scale by L 
and quantify it through: 



L - 



9e_ 
02 



(16) 



where H is the gas disc semi-thickness. It foflows that L is 
a small fraction of H. 

We now decompose the variables in equilibrium value 
and perturbation as follows: 



V = Vo + V 



AV = AVo + Av 



p = po{l + S) 
P = Po{x) + poh, 



(17) 
(18) 
(19) 
(20) 



where the zero subscripts refer to the unperturbed back- 
ground. Coordinates of the velocity perturbations are given 
by: V = uex + vBy + we^, Av = Aitea; + AvBy + Awe^. We 
then decompose the perturbation into axisymmetric Fourier 
modes oc exp [i{kxX + kzZ — ojt)]. Here, and kz are real 
wavenumbers and uj = ajsR -I- is is a complex frequency, with 
ujsn a wave frequency and s a growth rate (if s < 0, the 
perturbation is damped). 

YG05 studied the resulting sixth-order system numer- 
ically and found unconditional instability, regardless of the 
value of the dimensionless parameters Ts and fp, with a 
subdynamical growth rate scaling like ^Its in the tight- 
coupling limit. For a given Ts and fp, growth was max- 
imised in the k^-k^ plane in a long-wavelength ridge (with 
kz ~ fc^fstop/(i^/g)) and a short-wavelength branch with 
fcz ^ kx {kx having a preferred value). Growth was strongly 
suppressed in the former for dust-to-gas ratio near unity 
{fp ^ fg ^ 0.5). YG05 found that the wave speed uj^/k^ 
was bounded by [/sum defined as the sum of the background 
radial velocities of the dust and the gas, respectively, and 
growth rates (generally smaller than |fcj;|C/aum/2) were max- 
imised when the said wave speed was about [/aum/2. 



2.2 Terminal velocity approximation 

To gain physical insight on the streaming instability, 
we need some simplifications to make the problem tractable 
analytically. We shall adopt the "terminal velocity approxi- 
mation" (YG05), that is take the left-hand-side of equation 
1 10 1 to be zero, such that the relative drift is given by: 



AV : 



VP, 



^stop ■ 



(21) 



This approximation holds provided (i) the perturbation 
and dynamical timescales {lo~^ and are longer than 

fgtstop = PsCl/{pvt) (where the equality holds in the Ep- 
stein drag regime) and (ii) the lengthscale 1/fc of variation 
is longer than the "stopping length" /ggetstop = fgTsL (see 
Appendix B for justification). 

We will additionally neglect F and G because of the 
smallness of AV, an approximation YG05 find to be safe 
for Lj <^ Q,. 

The linearised system of equations is then: 



icuu — 2Qv + ge3 + ikxh — 



20. 



— ujw + kzh — 



ikxU + ikzW — iujS = 



(22) 



(23) 



(24) 



(25) 



ik^u + ikzW ~ ik^geifp - fg)tstopS + fpk^tstoph = (26) 

where zero subscripts have been dropped, k^ = k^ + k^, 
and we have introduced the epicyclic frequency k, given by 
= R-^d{R'^n^)/dR. In a Keplerian disc, k = Q, but 
one could also imagine studying streaming instabilities in 
any rotating fiuid, possibly even by direct experimentation. 
Equations (|22[)-([24| are derived from perturbation of the x, 
y and z components of equation ([9|. Equations ( 25 1 and ( 26 1 
stem from perturbing the continuity equations ([7|) and^l 
respectively. (In equation [26] we also make use of equation 
[21|.) Notice that the important term proportional to {fp — 



fg) in equation (26 1 derives from a sort of buoyancy force 



that arises in the two fiuid systems, with an effective density 
of f? I pp. When fp = fg, this effective density is stationary 
with respect to perturbations in pp. 

One can straightforwardly obtain a dispersion relation 
by setting the determinant of the above 5x5 system to zercj^ 



ifptstopi^* +0J^ + {ifpfi^ + kxgefg)tstopOj'' — ( K- 



~\~kx 



9et,top{fp ~ fg) ^ (27) 



Series solutions in Ts carried out by YG05 have shown that 
the roots fall in three branches, two epicycles giving rise to 
damping, and a secular mode, with: 



<^3i ~ kx{fp — fg)getstop + o(tstop) 



S = ifpt 



p'-stop 



(ifp - fg)ge^^ + o(tLp), 



(28) 



(29) 



^ We note that the first term of the quadratic coeflicient differs 
from that of equation (39) of YG05. While we believe that equa- 
tion \27\ of the current paper is correct (equation \29\ corrects 
their equation (44) accordingly), the YG05 error appears to be 
typographical and in no event does it afi'ect their exact results. 
Also, YG05 drop the quartic term in their equation (39); this is 
however consistent with cjfgtop 1- 



4 Jacquet et al. 



to leading order in the stopping tim^ (The notation o{X) 
indicates that the ratio o{X)/X tends to zero as X — )■ 0). 
It can be seen that feedback is essential to the instability, 
as ignoring it would amount to setting /p = everywhere. 
Also, the necessity of the background drift as measured by 
Qs is evident. More rigorously, it may be seen that for = 0, 
the dispersion relation becomes, if we disregard the u = 
mode: 



fplitstopX^ — X + fpKtstopX — 



= 



(30) 



with X = iuj/n. One can show that, for sufficiently small 
Ktstop, this cubic has three real positive roots (in terms of 
X), corresponding to lj being a purely imaginary number of 
negative imaginary part . 



2.3 The secular mode in a reduced system 

It is possible to simplify the system further while re- 
taining the leading-order physics of the instability. First, we 
note that, in general, uj <^ Q, as may be judged from the 



exact results of YG05 or equation ( 28 1 , for wavelengths that 



are not significantly shorter than the radial pressure scale L. 



Equation ( 23 1 gives us 



It is then easy to show that in equation ( 22 1 the first term 
(the acceleration term —iuju) can be neglected. It follows 
that, to leading order, the x and y equations of motion of 
the centre of mass relax to a form of geostrophic balance, 
i.e. Coriolis forces effectively cancel out the pressure gradi- 
ent and the buoyancy force. This is important, as it means 
that pressure perturbations do not effectively drive centre- 
of-mass motions in the orbital plane. So this crucial stabil- 
ising tendency is consigned to a subdominant role. It also 
means that the centre of mass executes modified epicycles 
in the orbital plane, to leading order. 

Second, we neglect the u perturbation in the two con- 



tinuity equations ( 25 1 and ( 26 1 which we are permitted to 



do if we restrict ourselves to long radial wavelengths. From 



equations (22l-(23l we can obtain the following scalings 



max(fc^^ h,uj^ &)■ 



Combining these with (|28fc and equation (241, we see that 
kxU is subdominant hereijj 

kx L ^ 1, and 



)( 



(31) 



which can be satisfied if k^H ~ 1 (see equation ( 16 1) and if 

ka:/k, < 1. 

^ It must be cautioned that the terminal velocity approximation 
discards terms of third and higher order in tstop that might con- 
tribute to the leading order expansion of the growth rate and 
hence could in principle affect its sign (which the YG05 calcula- 
tions show not to be the case here). Since we are here interested 
in understanding which ingredients give rise to instabiUty in a 
model system, rather than studying its exact properties (as in 
YG05), we content ourselves with this heuristic approach. 

The second condition results from requiring 
kzW with h eliminated from equation (|24|. 



These approximations give us a simplified set of dynam- 
ical equations: the equations for the vertical centre of mass 
velocity, particle density conservation, and the incompress- 
ibility of the gas, in which drag effects enter implicitly. These 
are now 



(32) 
(33) 
(34) 



— ifptstop yjr'^ j + + kxgetstopifg ~ fp) = 0, (35) 

from which we recover the secular mode and its growth rate 



— iujw + ikzh — 

ikzW — iioS — 

ikzW — ik^geifp - fg)tstopS + fpk^Utoph = 
yielding the following quadratic dispersion relation 
k 



w 




Lj/kz 


5 


= S 


1 


h 




_ {io/kzf _ 



as given in equation (29l! The corresponding eigenvector 
reads: 



(36) 



The second root corresponds to a damping s — 
— {k/kzY /(/pistop) but violates the condition a;4top «C 1 
of validity of the terminal velocity approximation and will 
thus not be discussed further here. 

It is also possible, if less rigorous, to see this reduc- 
tion dire ctly by inspection of the quartic dispersion rela- 

^ < 1 and fc^gTe/ft^ < 
one is allowed 



tion (271. If we have kj:ge.tl^op(kx / k 



1 (w 

to modify the quartic term as 



rich imply the conditions of (31 1), 

ifptstop {k/kz 



ifp (kx/kz) fstopW <^ u! and the quadratic one as {ifpK + 
kxgeifg - fp))tstop^^ since kxgefp < i/pK^. Under these 
conditions, the quartic may be factored as: 



'^^ + + kxgetstopifg - fp) 



ifptstop ( , 



such that we retrieve the quadratic dispersion relation of the 
reduced system. 

Hence, it appears that the mathematical essence of the 
instability can be isolated if we make the following set of 
assumptions: (a) the terminal velocity approximation, in 
which stopping time is short and so the relative velocity 
is determined from the steady balance (21 1; (b) geostrophic 
balance holds in the remaining momentum equations, which 
account for the horizontal centre of mass velocity (the cen- 
tre of mass executes modified epicycles in the orbital plane); 
and (c) that radial wavelengths of perturbed quantities are 
long (of order L or longer). In so doing, we have reduced the 
order of the system from 6 to 2. 

The ingredients for instability, apparently, are: a verti- 
cal velocity generated by a vertical pressure gradient (equa- 
tion [32]); accumulation of particles by the associated ver- 



tical flux of background particles (equation 33 ); and fi- 
nally gas incompressibility (equation [34| ). This last equa- 
tion, though difficult to interpret, is also the location that 



drag forces appear explicitly, via (211. Lastly, we emphasise 
again the importance of geostrophic balance, which offsets 
pressure gradients by Coriolis circulation, rather than radial 
flow. Thus radial pressure gradients do not drive stabilising 
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radial motions wiiich might alleviate the self-same gradients 
(but note that vertical pressure gradients remain in full). 
Therefore, rotation may not make an explicit appearance 
in these equations but its influence is crucial, as the next 
section demonstrates. 



3 STREAMING STABILITY IN THE 
ABSENCE OF ROTATION 

To emphasise the importance of rotation, and 
geostrophic balance in particular, let us consider the same 
problem with neither Coriolis force nor background shear, 
that is, in a nonrotating frame, a uniform flow of both dust 
and fluid, drifting relative to each other because of a pressure 
gradient, which is compensated for by an external gravita- 
tional field (and/or an inertial, non-Coriolis acceleration). 
It is actually more straightforward to handle this problem 
directly in terms of the dust and solid perturbations, with 
none of the approximations used in §3. We denote by Vp,g 
the perturbations of Vp,g, Sp = S/ fp the logarithmic density 
perturbation of the gas, and hg = h/fg. 

Perturbing equations ([l|-(|4]l (with = 0), still assum- 
ing incompressibility, gives the following system: 



{-i{uj - k ■ Vp) + 



stop 



-)Vp 



''Stop 



f-stop t-stop 



t'stop 



-Sp{Vp - Vg) + ihgk = 



(37) 



(38) 



ik ■ Vp — iujSp = 



(39) 
(40) 

where e = Pp/pg. identity matrix. Setting the determinant 
of the above 8x8 system to zero, the following dispersion 
relation results: 



a; - k • Vp + 



a;' + ( 



fgt stop 



k-(Vp + VJ)a; 



+ (k.Vp)(k.V,)-z 



. k- V 



(41) 



fgtstop , 

The first factor corresponds to a mode where the gas' ve- 
locity field is unperturbed, with the gas compensating the 
disturbance of the dust's with variation of pressure and den- 
sity. The quadratic factor (squared) has roots: 



i(k.(Vp + V,) + — ^ 

^ \ J 9 ''Stop 



1 ± 



V{1 + tie - l)(/,t,iop)2k ■ g,)2 - 4e(k ■ Seifgt.toprr^ ) , (42) 

with ge = VP/ p. Using the lemma | Re{^{l + ia)^ - h^) |< 
1 for (a,fe) G R2 (YG05), their imaginary part (s) is always 
negative. For k-ge(/gtstop)^ ^ 1 (as in the terminal velocity 
approximation limit), these roots may be expanded as 



k • (/pV, + /,Vp) - ie(k . ge)'(/,t.top)' (43) 



the latter corresponding to a damping where dust and fluid 
converge in concert toward equilibrium. 

There is thus no streaming instability in the absence 
of rotation (as was also clear from the YG05 numerical cal- 
culations). Rotation, as well as feedback, are both essential 
for the streaming instability, hence the difficulty in finding 
a simple toy model to explain it. 



4 TOWARD AN INTERPRETATION OF THE 
STREAMING INSTABILITY 

We now develop a framework with which to interpret 
the streaming instability. 

4.1 Particle concentration in pressure maxima 

We start with the equation governing the evolution of 
the dust mass fraction. Combining equations ([T]) and ([2| (we 
use here the full equations rather than those of the reduced 
system), we obtain: 

v-(p.(Vp-v .)) 
p p 

It is clear from this equation that dust is advected by Vp 
but is also subject to a mass flux associated with the rela- 
tive drift (the second term on the left-hand-side of equation 
( |45[ )). We might expect then to observe dust clumps corre- 
lated with AV (as noted by YG05 in the caption of their 
figure 6). In the terminal velocity approximation, equation 
(|45|) becomes: 



= 0. 



(45) 



Pp 



PgU 



' p J 



p 



AV • Vlnp - Utop 



VP 



(46) 



where the latter equality assumes the constancy of Pgtstop- 
There are two terms on the right-hand-side: one due to den- 
sity variation, and one to pressure variation. We shall show 
below that the latter is responsible for growth, not the for- 



mer (as was also ruled out by YJ07). Equation (46 1 may be 
rewritten in the useful form: 



^ + (/pV, + fgVp) ■ In^ 



Pgietop y2p 



(47) 



in which we have used incompressibility. Note that we have 
not yet linearised the system — this equation holds equally 
well in the nonlinear regime — and also that in this subsection 
we have made no hypothesis on the centre of mass dynam- 
ics (e.g. presence of Coriolis forces), since we have only used 
mass conservation equations and the terminal velocity ap- 
proximation. Its linearised form, however, with the Fourier 
dependence exp [i{kxX + kzZ ~ i^t)], yields 

P' 



k- (/pV, + /,Vp) + i/pfc^te 



p'p 



(48) 



with P' = ph and p'p = p5 the perturbation in pressure and 
(particle) density, respectively. 

What we have in (471 is a form of 'advective-reaction' 



and 



equation. Dust is advected at a velocity intermediate be- 



k- V 



fgts 



(44) 



/pV, + /,Vp = Vp + (V, - V) 



(49) 
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which corresponds to the C/sum of YG05, as in the frame 
used, k • V — 0. On the right hand side we have a source 
term proportional to the pressure Laplacian. This term will 
tend to draw dust towards pressure maxima in the x — z 
plane and away from pressure minima. Clearly, instability 
is intimately connected with this term, and the unstable 
mode works by concentrating dust at pressure maxima. We 
discuss both processes in turn. 



4.2 Advection term 

The left-hand side of equation ( |47[ ) gives rise to the 
streaming instability wave character: a density pattern in x 
and z will be advected at the velocity 



V, + (V, - V) 



(50) 



and in the linear regime this will be entirely radial. The ad- 
vection velocity is the sum of two parts: the particle velocity 
and a secondary contribution issuing from the relative drag 
(the bracketed term). The first term should be familiar, and 
simply describes the Lagrangian advection of particles by 
the particle velocity itself. The second term is novel. Phys- 
ically, it represents the fact that particle density maxima 
induce a decrease in the relative velocity - because the re- 
ciprocal drag of the two fluids is greater at that location 
(cf. the first term in equation 
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This decrease leads to 
an additional mass flux, which is anticorrelated with pp. An 
outward effective advection is the result, which is in addi- 
tion (and opposition) to that of Vp. Thus a density peak 
will tend to be both (a) pushed inward simply by the pri- 
mary particle velocity Vp, and (b) pushed outward by the 
second drag-induced mass flux that it has itself excited. It 
is easy to show that for fg = fp the sum of the two ve- 
locities is precisely zero, i.e. the two drift velocities cancel 
perfectly. Conversely, when there are very few particles, i.e. 
/p — ^ 0, the advection velocity is simply Vp and the pattern 
is carried by the particle fluid alone. This is because, in this 
limit, the particles' influence on the gas is so minor that Vg 
is efi'ectively zero. On the other hand, when there is virtu- 
ally no gas, fg — >■ 0, the advection approaches zero, because 
Vg — >■ (as well as Vp) in this limit. 



4.3 Instability term 

We have already mentioned that instability was related 
to pressure maxima (or technically, pressure perturbation 
maxima, as there is a background gradient, which however 
has zero Laplacian ). This is reminiscent of the concentration 
mechanism of test particles in long-lived, two-dimensional 
vortices (e.g. Barge and Sommeria|1995 Klahr and Boden 



heimer|2006 1, which in a Keplerian flow are anticyclonic 

In the streaming instability, this concentration is an ac- 
tive process, in the sense that feedback participates in the 
preservation of the high-pressure, particle-trapping zones. 
How does the instability grow? Schematically, as we have 
said, pressure maxima attract dust particles, as can be seen 
from the terminal velocity approximation given by equation 
(211. Then, dust drags the gas toward these same maxima 
(see the gas Euler equation Q, where the terminal velocity 
approximation ( |21[ | may be injected). This creates conver- 
gent flows that strengthen these pressure maxima. These in 




' Dust grains 

PRESSURE MAXIMUM ATTRACTS DUST ^ DUST DRAGS THE GAS PRESSURE MAXIMUM STRENGTHENED 

Figure 1. Cartoon of the proposed interpretation of the stream- 
ing instability (seen along the radial axis). Dust is drifting ra- 
dially inward. The pressure profile is drawn in blue (we ignore 
the background gradient in the picture) . We start with a pressure 
perturbation maximum. This attracts dust, which in turn drags 
the gas, creating convergent flows that strengthen the pressure 
maximum. Note the role of geostrophic balance in maintaining a 
circulation (schematically shown in magenta) which supports the 
pressure maximum. 



turn can attract more dust particles, and the process thus 
runs away. This is sketched in Figurell] 

Indeed, inspection of equation ( |48[ ) shows that growth 
is correlated with fpP'/p'p- /p is a measure of the feed- 
back of dust on the gas, while P' / p'p measures how well 
dust density and pressure maxima correlate with each other, 
or in other words, the efficiency of the dust-driven pres- 
sure loading. Growth is positive if both pressure and den- 
sity perturbations are not out of phase by more than 
90°. For example, YJ07's eigenvector "linA" has a pres- 
sure phase relative to particle density of 43 ° and a com- 
plex frequency uj = (-0.3480127 -f 0.4190204i)f2 while that 
same phase for "linB" is — 95 ° with a complex frequency 
uj = (0.4998786 + 0.0145764i)nQ For the secular mode in 
the reduced model, P' I p'p ~ h/5 — {ijj/kz)^ (see equation 
(36l) is real and positive to leading order. 

While the above identifies the positive feedback loop 
leading to instability, we have yet to account for two other 
physical ingredients that are necessary condition to it. In- 
deed, drag of the dust on the gas is actually insufficient alone 
to offset the pressure force that tends to accelerate gas away 
from pressure maxima. In the terminal velocity approxima- 
tion, we have: 



Pp Vp 



P9 



VP 



VP 

P 



(51) 



The first physical ingredient we wish to mention is rota- 
tional dynamics, whose necessity has been shown in §3. We 
have already mentioned in §2.3 that a geostrophic balance is 
established, by which Coriolis forces essentially balance the 
pressure force. The other physical ingredient is the back- 
ground solid-to-gas drift (due to the pressure gradient), as 
noted in §2.2. Mathematically, a nonzero drift is required 
for particle density to affect the dynamics, as, in the re- 
duced system, this can only happen through the "buoyancy 
force term" (proportional to the pressure gradient) in the 
gas continuity equation (|8| as mentioned in §2.2, which con- 
tributes to the advection term discussed in §4.2. We suggest 



* s is slightly positive, instead of slightly negative, presumably 
because of finite gas compressibility and/or corrections to the 



terminal velocity approximation. 
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that the physical explanation for this is that relative drift 
between the pressure maximum and the dust allows the for- 
mer to sweep a (radial) "headwind" of particles and be an 
effective "trap" for the dust. 

We note that the processes discussed in the previ- 
ous paragraph are restricted to the x-y plane and empha- 
sise radial drifts and gradients. The fact that the reduced 
model seems to emphasise vertical gradients — although the 



leading-order growth rate given by equation ( 29 1 is nega- 
tively correlated with — may be interpreted as resulting 
from incompressibility, as radial velocity perturbations of 
the gas must then give rise to vertical perturbations and 
gradients (see e.g. flow pattern in Fig. 5 of YG05). 



4.4 Streaming instability as collective drag? 

The above demonstration leads us to question the cor- 
respondence suggested by CYIO between the streaming in- 
stability and the toy model of [Goodman and Pindor| ( |2000[ ), 
which was originally designed for a dust layer-scale insta- 
bility. In this one-dimensional, single-fluid model, the ac- 
celeration is the sum of g, a proxy for gravity, pressure and 
inertial forces, and —i/d{'S)v a density-dependent drag accel- 
eration proportional to the velocity v of the fluid. As long 
as dvd/dT, ^ 0, linear analysis shows the system to be over- 
stable (CYIO). 

The fact that this toy model is one-dimensional ap- 
pears to prevent it from capturing the specificities of Corio- 
lis forces, unlike other gravitational or inertial terms. It may 
also be questionable whether in the linear phase, and in the 
local treatment of the instability, the drag can be considered 
to be collective in the sense of the above toy model, and 
what would play the role of the (single) fluid of this model. 
Certainly, the latter cannot be the gas+dust fluid, as drag 
is internal to it (see equation |9]). But even if we take the 
dust for example, and use the terminal velocity approxima- 
tion (or higher-order corrections) as a proxy to eliminate the 
gas velocity (we cannot treat it as tmposed — this is indeed 
the point of the streaming instability), the drag accelera- 
tion becomes —VP/p which makes no reference to the dust 
velocity. Same holds for the gas. Finally, as we have noted 
earlier, the variation of the relative drift with density is not 
responsible for growth in this approximation and hence does 
not conform to the notion of a collective drag underlying the 
instability. 



5 CONCLUSION 

Through successive approximations, we have obtained 
a reduced system giving rise to the streaming instability. 
It would appear that the essence of the instability can be 
effectively established given: 

• the terminal velocity approximation 

• that a form of geostrophic balance holds for the planar 
centre of mass velocity 

• that radial wavelengths are sufficiently long. 

The leading order (in the stopping time) expansion of the 
growth rate is the same as the fourth-order system with the 
terminal velocity approximation. 



We interpret the streaming instability as arising from 
dust pileup in pressure perturbation maxima, with the dust 
then dragging the gas and hence strengthening the maxima. 
For this enhancement process to occur despite the pressure 
force that tends to "smear" these maxima out, it is necessary 
that there be a background solid-gas drift, and also Coriolis 
forces that help establish a geostrophic balance with the 
pressure force. The role of Coriolis forces would appear to 
indicate that one-dimensional toy models do not adequately 
capture the instability. 

Larger amplitude vertical flows will (when outside of 
the linear regime) transport density away from maxima and 
hence limit the growth of the instability, and would compete 
with other nonlinear effects like particle clumping ('traffic 
jams'). Incidentally, another process in the nonlinear satu- 
ration that will work against particle clumping is particle 
pressure, not usually modelled. The large local densities of 
nonlinear clumps will lead to enhanced collision frequencies 
and consequently to a particle pressure which will resist fur- 
ther concentration. Attempts to establish the conditions for 
self-gravitational collapse of such clumps should in principle 
include this effect. 
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APPENDIX A: APPROXIMATING DUST AS A 
PRESSURELESS FLUID COUPLED TO THE 
GAS 

It is of interest to investigate in some generality how 
the description of dust as a pressureless fluid, used in many 
simulations, arises. As CYIO point out, the usual criteria 
for gas molecules do not apply. Intuitively, one expects a 
fluid description to be viable if the stopping time is short 



in some sense (YG05, see also YJ07). Here, we verify this 



intuition using a kinetic approach. Garaud et al. ( 2004 1 ad- 



dressed this problem for one-dimensional settling in a static, 
homogeneous gas, finding that velocity dispersion would be 
quickly damped for tightly coupled particles in the Epstein 
regime. 

We introduce a distribution function /(r,Vp) of the 
number of particles in the (r-Vp) phase space. The conti- 
nuity equation in this space may be written as: 



dl 
dt 



+ Vp 



dr dvr, 



(Al) 



with Fp the total force exerted on one particle (except 
particle-particle interactions), mp the particle mass and 7coii 
the collision integral. We have used • (/vp) = Vp • ^ 
but here, • Fp 7^ because Fp includes the drag force 
—rrip (vp — Vg) /istop, hence the departure of the left-hand- 
side from the standard Boltzmann equation form[^ 

Provided collisions conserve the total number and mo- 
mentum of particles, integration over velocities of equation 
( Al I mutliplied by Vp, yields: 

1 



VpVp = 



Pv 



(A2) 



with Vp = (vp) the mean particle velocity, S = 
/9p((vp — Vp) (vp — Vp)) the stress tensor. For any particle- 
wise quantity x, we have defined: 



{x) 



Pv 



xfd^Vp 



(A3) 



Treating the dust as a pressureless fluid amounts to neglect- 



ing the stress tensor in equation ( A2 1. To evaluate this, it is 



of interest to derive its evolution equation: 



rrip / /coii(fp - V^p)(wp - V/)d^Vp 



- 1 [f;{vI - v^) + F^(v; - v;)) fd\p 



(vp - Vp) 



or 



v;){vi,-v^)d\p. 



(A4) 



The first term on the right-hand-side relates to the redis- 
tribution and loss of relative velocities during collisions, the 
second to the action of velocity-dependent forces, and the 
third to transport of velocity dispersion between neighbour- 
ing "fiuid elements". The drag contribution to the second 
term is — 2S'-'/4top (for a velocity-independent tstop)- 

In the absence of particle-particle interaction (/coii = 0), 



integration over space of equation ( A4 \ yields 



+6jMn'{5")v + {|^5'=')v)(A5) 



{.■■)v denotes a spatial averaging, eijk is the Levi-Civita ten- 
sor and ft is the instantaneous rotation vector of the refer- 
ence frame. Damping by gas drag dominates the evolution 
of the velocity dispersion if tstop is shorter than f2~^ and the 



^ Thus, technically, what jCaraud et aZ. | ( |2004^ called an interac- 
tion term, although their approach was collision-free, actually is 
the corresponding correction to the left-hand side. 
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viscous heating rate |VVp| ^. This is thus the criterion for 
validity of the pressureless fluid approximation. 



APPENDIX B: VALIDITY OF THE TERMINAL 
VELOCITY APPROXIMATION 

Perturbing equation |To| yields: 



— icjAv — 2nAvex H Auey 

2 

-ik^Qetstop (v + {fg - fp)Av - fgAYS) 

Av + SAY ft , , 
= +i^k (Bl) 

,/gIstop Jg 

(This can also be obtained by injecting equation 31 into 
equation 27 of YG05). The terminal velocity approximation 
amounts to setting the left-hand-side equal to zero. 

The constraint ujfgtstop ^ 1 results from comparing 
the first term on the left-hand-side with the first one on the 
right-hand-side. The latter dominates the next two terms 
on the left-hand-side if fi/gtstop ^ 1, and the fourth and 
fifth ones if /gk^Qetstop ^ 1 (recall that both fp and fg are 
smaller than 1). To justify the assertion for the fourth term, 
we combine the perturbations of ([7| and ([sjl (equations 28 
and 29 of YG05) to obtain: 

k-v^ ^f,^: (B2) 

and use the expectation that the wave speed is of order the 
background drift velocity, as verified by YG05|^ The final 
term on the left-hand-side is dominated by the second one 
on the right-hand-side under the same condition. 

The condition /gtstop is implied anyway 

for the validity of the fluid approximation (see Appendix 
A). Since uj < Q, the condition on Ts (^ l//g) is the most 
stringent one. These conditions are somewhat less stringent 
than those mentioned by YG05 (which are thus sufficient), 
and more symmetric with respect to dust and gas (as they 
converge in concert toward the terminal (relative drift) ve- 
locity because of their mutual interaction). 



° Actually, equation |[B2| does not really con strain v, but if we 
adopt the scaling v ~ from equation |23| , the condition for 
this component reduces to the already obtamcd flfgistop <C 1. 



